On the Profiles and the Polarization of 
Raman Scattered Emission Lines in Symbiotic Stars 
II. Numerical Simulations 



K. W. Lee^ & Hee-Won Lee^ 



^Dept. of Astronomy and Atmospheric Sciences 
Kyungpook National University, Taegu, Korea 

^ Research Institute for Basic Sciences and Dept. of Astronomy 
Seoul National University, Seoul, Korea 

email: hwlee@astro.snu.ac.kr 

ABSTRACT 

A Monte Carlo method is used to calculate the profiles and the polarization of the 
Raman scattered O VI lines (AA6827, 7088) in symbiotic stars, which are believed to be 
a binary system of a cool giant and a hot star with an emission nebula around it. A 
point-like isotropic UV radiation source is assumed and a simple spherical wind model is 
adopted for the kinematics of the scattering material from the cool giant. 

We first investigate the case where the incident line photons are described by a Gaus- 
sian profile having a width of 10^ K. We subsequently investigate the effects of the ex- 
tended ionized region and non-spherical wind models including a disk-type wind and a 
bipolar wind. The cases where the emission source is described by non-Gaussian profiles 
are briefiy studied. 

Finally as an additional component for the kinematics of symbiotic stars the orbital 
motion of the hot component around the cool giant is included and the effect on the 
spectropolarimetry is investigated. In this case the polarization direction changes around 
the red part of the Raman-scattered emission lines, when the observer's line of sight is 
perpendicular to the orbital plane, and no such effect is seen when the line of sight lies 
in the orbital plane. Furthermore, complex peak structures are seen in the degree of 
polarization and the polarized fiux, which have often been observed in several symbiotic 
systems including RR Tel. 

Brief observational consequences and predictions are discussed in relation to the 
present and future spectropolarimetry for symbiotic stars. It is concluded that spectropo- 
larimetry may provide a powerful diagnostic of the physical conditions of symbiotic stars. 



1. Introduction 



A symbiotic star is believed to be a binary system consisting of a cool giant and a 
hot star providing H-ionizing photons for an emission nebula around it (Iben & Tutukov 
1996). It is expected that the hot component is characterized by a typical temperature 
T ~ 10^~^ K and that the binary motion has a typical period of P ~ 10^"'^ days. 
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The observed spectral parts of the emission in the symbiotic systems cover a very 
large range including radio, IR, optical, UV and X-rays. The complicated profiles and the 
significant variabilities often seen in the UV emission lines imply that the dynamics and the 
physical conditions are not simple enough to be described by a unique model {e.g. Vogel 
Sz Nussbaumer 1994). 

It has been known that a significant fraction of symbiotic stars exhibit broad emis- 
sion features around A6827 and A7088 with widths up to 20 A, which are an order of 
magnitude larger than those of other emission lines in their typical spectra {e.g. Allen 
1980). These broad features are identified by H. Schmid (1989) as the Raman scattered 
O VI AA1032, 1038 by hydrogen atoms. During the scattering process, the atom in the 
ground state absorbs the incident line photon with frequency Ui to be excited to an inter- 
mediate state followed by de-excitation to the 2s state emitting a photon with frequency 
Vf — i>i — ff^-i where fa is the frequency of the Ly a transition. 

The Raman scattering nature explains the large width of the features and many ob- 
servational characteristics of symbiotic stars such as the existence of highly ionized lines. 
Observational confirmations of the Raman scattering mechanism include the spectropo- 
larimetry which shows that the features are highly polarized in contrast with other negligi- 
bly polarized lines {e.g. Schmid & Schild 1994). The identification of the Raman scattering 
nature is corroborated by the UV observation by Espey et al. (1995), who showed that the 
symbiotic star RR Tel exhibits a very strong O VI doublet with the intense broad features 
around A6827 and A7088. 

Valuable kinematic information on the emission nebula of a symbiotic system is ob- 
tained in a detailed analysis of the profiles of the emission lines, which often accompany 
non-trivial structures including double-peaked or partial absorption troughs {e.g. Mueller 
&; Nussbaumer 1985, Pereira et al. 1995). The Raman scattered features can provide a 
useful diagnostic to constrain many kinematical parameters of a symbiotic system, because 
they have broadened profiles due to the relative motion of the stellar wind around the cool 
giant with respect to the emission nebula. Furthermore, the width of the profiles is also 
enhanced by the inelastic nature of the scattering by almost an order of magnitude. 

Recently Harries & Howarth (1996, 1997) provided a good amount of spectropolari- 
metric data and performed a fairly complete Monte Carlo calculation about the Raman 
scattered emission features (see also Schmid 1996). In their observational data some sym- 
biotic stars show the polarization direction fiip around the red part of the Raman scattered 
feature. This kind of behavior usually accompanies multiple peak structures in the polar- 
ized flux {e.g. Espey et al. 1995, Schmid 1996). However, in the spherical wind model, 
many Monte Carlo computations show that the polarization flip usually occurs around the 
center of the feature calling a need to investigate non-spherical models. 

This is the second paper in a series on the polarization and the line profiles of the 
Raman scattered flux. In Lee & Lee 1997 (hereafter Paper I), we reviewed the basic atomic 
physics of the Raman scattering of incident photons shortward of Ly a by hydrogen atoms, 
and computed the profiles and the polarization for a few cases using the single scattering 
approximation (see also Lee & Lee 1996). A more complete review about the atomic 
physics is provided by Isliker et al. (1989) and the references therein. In this paper we 
calculate the line profiles and the polarization of the Raman scattered O VI lines for various 
kinematic models using a Monte Carlo method. 



2 



The paper is composed as follows. In section 2, the basic procedures of the Monte 
Carlo method and the model descriptions are given. Using the Monte Carlo code we 
investigate some fundamental properties of the Rayleigh-Raman scattering processes and 
discuss them in section 3. The results of the numerical calculations are presented in section 
4. In section 5, observational implications are commented. Finally, in section 6, a summary 
is presented with discussions about the significance of spectropolarimetry in understanding 
the symbiotic systems. 

2. Model 

There has been much observational work on symbiotic systems, showing that compli- 
cated kinematics is responsible for variabilities and complex structures in the line profiles 
{e.g. Mueller & Nussbaumer 1985, Pereira et al. 1995). We do not attempt to give a 
detailed dynamical model in this work, but adopt rather simple models in order to find 
out the main features contributing to the polarization structures and the profiles in the 
Raman scattered lines {e.g. Harries &: Howarth 1996, Schmid 1996). In particular, the 
polarization behavior is expected to depend upon the scattering optical depth of the sys- 
tem, which naturally measures the number of scatterings before a given photon gets into 
the line of sight. The Raman scattering of the O VI photons is also characterized by an 
enhanced Doppler effect, according to which the wavelength shift of an incident photon 
increases by a factor of ~ 7. Therefore, the kinematics of the scatterers may combine with 
the scattering optical depth structure to contribute to the complicated behavior of the 
polarization of the Raman scattered features. 

In this section, we summarize the basic points of the Monte Carlo code and the 
kinematic models adopted in our numerical simulations. For the Raman and Rayleigh 
scattering cross sections of the O VI doublets, we use the results in Paper I, in which we 
showed that 

(7fla^(1032) = 7.5 ar , cTi?am(1038) = 2.5 ar 

CTi?aj; (1032) = 34 CTT , CTjjaj, (1038) = 6.8 (It 

where ctt = 6.6 x 10~^^ cm^ is the Thomson scattering cross section. 

2.1 Monte Carlo Method 

In this subsection we briefly describe the Monte Carlo method to compute the profiles 
and the polarization of the Raman scattered O VI doublet lines. It is first assumed that 
there is an emission region embedded in the scattering medium consisting of the extended 
atmosphere of the cool giant. The incident UV line photons are further assumed to be 
emitted isotropically and unpolarized. We start with the incident photons governed by 
a Gaussian profile corresponding to a Maxwell-Boltzmann distribution. We subsequently 
investigate the cases where the incident radiation field is described by synthetic profiles 
such as double-peaked profiles often seen in the observational data. 

In Fig. 1 we show a schematic geometry of a symbiotic star system, which is adopted 
for our Monte Carlo calculation. The coordinate system is chosen in such a way that the 
binary axis coincides with the 2;-axis and the observer's line of sight is the a;-axis. The hot 
star is represented by the circle on the right side, and the bigger circle in the origin is the 
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cool giant. A spherically symmetric stellar wind is depicted by radial arrows around the 
cool giant. 

For a numerical simulation of the Rayleigli- Raman scattering it is essential to compute 
the free path length at a position for a given wave vector and wavelength. The mean free 
path of a UV line photon is given 

If = {natoty\ (2.2) 
where the total scattering cross section atot is defined by 

O'tot = CTRam + Ray (2.3) 

A given photon travels up to a distance I with a probability 

Pi = l- exp{-l/lf) = 1 - exp(-TO, (2.4) 

where 

Ti = l/lf = j ds n{s)(Ttot- (2.5) 

In terms of a random uniform variable r in an interval between and 1 , t/ can be generated 
by a transformation 

T« = -ln(l-r), (2.6) 

identifying r with pi. The inverse transformation of Eq. (2.5) gives the free path length of 
the incident photon. 

A photon travels a free path I before a hydrogen atom Raman-scatters it with a 
branching ratio rRam = o'Ram/f^tot- If a random number is greater than rRarm then we 
regard this scattering as Rayleigh, and otherwise the photon is Raman-scattered. 

When the scattering type is Rayleigh, the escape condition should be examined before 
the next scattering site is determined. The escape condition is simply 

n > Tesc, (2.7) 

where Tesc is the scattering optical depth to an observer at infinity, i.e. , 

poo 

Tesc = ds n(s)atot- (2.8) 



If the escape condition is met, then the UV photon is regarded as being emitted to infinity. 
Otherwise a new wave vector and the polarization components are calculated to be used 
for the subsequent scattering. 

On the other hand, if the photon is Raman-scattered, then we record the photon 
with the polarization according to the wavelength. Here, we assume that the continuum 
absorption around Ha is negligible and that the Raman-scattered photon escapes freely. 
The effect of continuum opacity has been extensively investigated by Schmid (1992, 1996) 
and Harries Sz Howarth (1997). In order to keep the standard deviation no larger than 1.5 
percent, we generate typically 5 x 10^ photons for each run. 
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In the following section, we discuss the stellar wind models and calculate the relevant 
scattering optical depths and the escape conditions. 

2.2 Stellar Wind Models 

We discussed the kinematic models in Paper I, which we adopt and again introduce 
in this paper (see also Nussbaumer & Vogel 1987). A typical wind velocity law can be 
written as 

v(r)=^oo(l-i?./rf f 

where Voo is the terminal wind velocity, i?* is the radius of the cool giant, p = r / R^, is 
the radial distance in the unit of i?*, and /3 is a positive constant. The density profile 
corresponding to the velocity law Eq. (2.9) is given by 

n(r) = nop-'(l-p-^)-^, (2.10) 

where no is a typical density given in terms of the mass loss rate M, and the proton mass 
rup by 

no = M /AirRlmpVoo- (2.11) 

We again choose /? = 1 as in Paper I. 

Therefore, the total scattering optical depth between position 1 and position 2 is given 

by 

ds n(r) atot 



= To 



1 y^^^(p-l) 



where s is the running parameter along the photon path, and the representative scattering 
optical depth 

To = noR^atot- (2.13) 

Here, b is the impact parameter of the photon path with respect to the center of the cool 
giant divided by i?* (see Fig. 1 for a schematic geometry). In this paper, it is understood 
that all the distances are measured in units of the radius of the cool giant i?* . 

If the i-th scattering site is given by pi and the wave vector by k^, the next scattering 
site pi+i characterized by the scattering optical depth Tj from pi is obviously given by 

Pi+i = Pi + \pi+i - Pi\K ^2 

= pi + Si{Ti)ki, 

where Si{Ti) is the distance between the two scattering sites, which is dependent on the 
vectors pi, k^ and Tj. 

The functional form of Si can be written separately for the case where the impact 
parameter of the photon trajectory is larger than the radius of the cool giant and for the 
case where it is smaller than the radius of the giant, i.e. , either 6 > 1 or 6 < 1. 
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In order to express the functional form of Si(Tj), we perform an integration of 
Eq. (2.12). Introducing a function f{p) defined by 



m 

we have an equivalent relation 

1 



which gives p after /(p) is obtained. 
If 6 < 1, then we obtain 



f{pi+i) = coth 



coth ^ f{pi) - 



2 To 



for • r > 0, and 



coth 



coth V(Pi) + 



2 To 



for f < 0. 

On the other hand, if 6 > 1, then the result is 



f{pi+i) = tan 



tan V(Pi) + 



To 



Si{ri) = ^pl^-b^-VW^ 



(2.15) 



1 + VW^lfip) + b'{{l + VW^\f{p)}-'] , (2.16) 



(2.17) 



(2.18) 



(2.19) 



for kj ■ f > 0. If 6 > 1 and k^ ■ f < 0, then we introduce the scattering optical depth T5 to 
the impact point, which is given by 



n = [tan-i /(pO - tan-^ /(6)] . 



(2.20) 



Then the final result is written as 

/(Pi+i) = tan 



tan ^ /(pi) - 



2 To 



(2.21) 



^Si{T,) = M^-^p\^^-h'^ 



for Tj < Tfe, and 



/(Pi+i) = tan 



tan"^ f{b) + 



2 TO 



(2.22) 



for Ti> n. 

The escape conditions are met when the scattering optical depth Ti > Tesc- Here, the 
escape scattering optical depth t^sc is obtained from Eqs. (2.17)-(2.22) by letting Si — > oo, 
i.e. , 

2to 



2ro 
2ro 



coth V(Pi), 



for 6 < 1, kj • r > 0; 



TT 



--tanV(P*)|, for 6 > 1, kj -f > 0; 



(2.23) 



TT 



-tan-V(&) , for 6> 1, ki -f < 0. 



No escape is possible for the case where b < 1 and k^ ■ f < 0, because the total optical 
depth to the surface of the cool giant diverges for all pi. This means that the photons are 
scattered off before they reach the surface of the cool giant. 

The final procedure of the Monte Carlo simulation is to determine the Doppler factors 
and the polarization associated with the scattered photon. The polarization state of an 
ensemble of photons is described by a density matrix {e.g. Berestetskii et al. 1971, Lee et 
al. 1994), where the Stokes parameters Q, U, and V correspond to the difference of the 
main diagonal elements, the real and the imaginary parts of the off-diagonal elements of 
the density matrix, respectively. The Rayleigh-Raman scattering is characterized by the 
Rayleigh phase function, and the wave vector and the polarization of the scattered photon 
are computed by the same way as in the case of the Thomson scattering. 

For the Raman scattering of the O VI doublets, the total wavelength shift is given by 



AXf _ A/ k, • v(r) ^ kf ■ v(r) 

A/ Aj c c 



(2.24) 



where k^ is the wave vector of the incident photon and k/ is that of the outgoing photon. 
Eq. (2.24) is also valid for Rayleigh scattering, where A/ = A^. It is clear that in the case 
of Raman scattering the Doppler factor becomes large by a factor of about 7 due to the 
motion of the scatterers relative to the emission source. We show the contours of constant 
total scattering optical depth and Doppler factors in Fig. 2, which were introduced and 
discussed in Paper I. 

3. Fundamental Properties of Raman Scattering 



In this section we discuss the characteristic features of the Raman scattering process 
using the Monte Carlo method. Schmid (1992, 1995, 1996) also discussed the basic prop- 
erties of the Rayleigh-Raman scattering process (see also Harries & Howarth 1997) and we 
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elaborate further the properties of the scattering process in a similar way to deal with the 
random walk process. 

3.1 Single Scattering Approximation 

In Paper I, we briefly mentioned the single scattering approximation, which applies to 
the case where the scattering region is characterized by a small scattering optical thickness. 
This formalism can also be useful as a check of our Monte Carlo code, in the sense that it 
is the limiting case when the total scattering optical depth tends to zero. 

In Fig. 3 we present a result from our Monte Carlo code for a scattering optical 
depth To = 0.5 and compare it with the corresponding result from the single scattering 
approximation. The adopted parameters are described in Paper I. The Monte Carlo result 
is in good agreement with the result from the single scattering approximation within one 
standard deviation shown by the error bars in the degree of polarization. 

In the following subsection we quantify the single scattering approximation and inves- 
tigate some fundamental properties of the Rayleigh-Raman scattering process. 

3.2 Reflection from a Slab 

A UV photon with a wavelength shorter than 1216 A can be scattered either by 
a Rayleigh process or by a Raman process from a hydrogen atom. Therefore, the UV 
photon may be Rayleigh-scattered several times before it is Raman-scattered. Under the 
assumption that a Raman-scattered photon has a very small optical depth, it will escape 
freely from the scattering region and may reach the observer. The ratio of the number of 
the emergent photons which are Raman-scattered to that of the total incident photons is 
expected to be dependent on the ratio of the Raman scattering cross section to the total 
scattering cross section and the scattering geometry such as the total optical depth of the 
medium. 

In this subsection, we discuss briefly some of the fundamental properties of the 
Rayleigh-Raman scattering process of UV line photons in a simple scattering medium. 

H. Schmid, who proposed the Raman process in symbiotic systems and performed 
pioneering works in this subject, presented the basic results in his papers {e.g. Schmid 
1992, Schmid 1996). A Monte Carlo simulation is particularly useful in describing the 
scattering process, because we can conceptually divide the emergent photons according to 
their scattering numbers and perform a detailed analysis. Pursuing in this line of reasoning, 
we collect the photons reflected from the both sides of a slab of flnite scattering optical 
depth, using the Monte Carlo code. It is assumed that the slab is illuminated from outside 
and that the distance from the slab to the incident raditation source is much larger than 
the size of the slab, so that the incident photons enter the slab effectively normally. The 
emergent photons are subsequently divided into their number of scatterings, and we record 
their flux and the degree of polarization for further analysis. 

Fig. 4 illustrates the result from the Monte Carlo calculation on the Rayleigh and 
Raman reflected components from a finite slab as a function of the number of scatterings. 
On the vertical axis is shown the logarithm of the number of the photons reflected from the 
slab. The total scattering optical depth Tg of the slab is chosen to be Tg = 0.5, 1, 5, 10. 
The ratio rjiay = crRay/(^tot of the Rayleigh scattering cross section to the total scattering 
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cross section is taken to be ruay = 0.2, 0.8, where the case rjiay = 0.2 is represented by 
the hght Unes and the thick Unas are for the r^ay = 0.8 case. 

When the scattering optical depth is small and rfiay = 0.8, the Rayleigh-scattered 
flux is larger than the Raman-scattered flux, whereas the converse is true for large scat- 
tering optical depths. This is because the Rayleigh-scattered photons are trapped as the 
scattering optical depth increases, whereas all the Raman-scattered photons are assumed 
to escape on the spot. 

As is usually expected the main contribution to the scattered flux is due to singly 
scattered photons. Both the flux and the degree of polarization show exponential decrease 
as a function of the number of scatterings. This behavior is the foundation of the single 
scattering approximation discussed in the previous section and in Paper I. 

We may give a semi-quantitative argument about the fraction of the Raman scattered 
flux to the total incident fliix. For a given slab of total scattering optical depth r^, a 
fraction (1 — e~'^'') will be scattered at least once. Let's denote the total incident number 
flux by /o. Then = /o(l — e""^") is the total scattered (number) flux. Let f{n) be the 
number flux scattered no less than n. From this deflnition, it is obvious that 

/(I) = (3.1) 

Furthermore, a fraction (l—ri^ay) of fs will be Raman scattered and escape the region. 
Hence if we denote by fRamip) the Raman scattered (number) fliix scattered only n times, 
then 

fRamil) = (1 - rRay)fs = (1 - rRay)f{\). (3.2) 

In fact, this relation holds for any n, that is. 



fRam{n) = (1 - rRay)f{n), (3.3) 

because any Raman scattered photon does not suffer a subsequent scattering. 

In a similar way we denote by fRayin) the Rayleigh scattered emergent number flux 
scattered only n times. Then we may write 



fRay{n) ^ rRay(3{n)f{n). (3.4) 

Here, the function /3(n) is loosely deflned as an escape probability from the n— th scattering 
site. Because the scattering process is similar to the random walk process, there is no 
deflnite site of n— th scattering, and we can only state in a probabilistic way. If we borrow 
from the solution of the random walk problem, we may state that the photon diffuses up 
to the position where the scattering optical depth from the surface is approximately y/n. 
Therefore, we may tentatively set 

P{n) = \{e-^ + e"^^^-^)). (3.5) 

When the number of scattering exceeds r|, then the above formula does not give any 
meaningful results. However, the flux /(n) decreases approximately exponentially, the 
result will not be sensitive to the choice of /?(n) for sufficiently large n. 
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Having defined fRamin), fRayin), and P{n), we can formally give a recursion formula 
for f{n), which is 

/(n + 1) = /(n) - fRamin) - fRay{n) 

= rRay[l-(3in)]fin). ^ ' ^ 

The functional dependence of (3{n) on the scattering optical depth and the scattering 
geometry determines the basic properties of the Rayleigh-Raman scattering process. 

In Fig. 5 we plot fRami'n), fRayin) as a function of the scattering number n, which 
are obtained recursively using Eq. (3.6). A comparison is made with the corresponding 
Monte Carlo results. The solid lines show the Monte Carlo results and the dotted lines give 
fRam{n) and fRay{n) obtained analytically. Eq. (3.5) is used for the escape probability 
for the upper panel. In the bottom panel we show a different set of fRamij^) and fRayin) 
using another formula for the escape probability, defined by 

/3'(n) = ^exp(-n«-^). (3.7) 

The adopted parameters are Tg = 10, and VRay — 0.8, which is relevant for the O VI 
doublets AA1032, 1038. Due to the choice of the parameters we see that the Raman- 
scattered fluxes are larger than the Rayleigh-scattcred coimterpart. As is shown in Fig. 5a, 
the agreement is very good except when the number of scatterings becomes large, where 
the small remainder flux does not give much signiflcance. Far better agreement is obtained 
in Fig. 5b. Considering the agreement shown in Fig. 5, the preceding argument may be 
said to closely describe the scattering process in the scattering medium. 
The total Raman scattered flux FRam is given by 

FRam = XI fRam{n) = (1 - TRay) ^ f{n). (3.8) 

n n 

Similarly, we may write the total Rayleigh scattered flux FRay 

FRay = ^fRay{n) = T Ray^ (3{n) f {u) . (3.9) 
n n 

Hence, we can estimate the ratio RRam of the total Raman scattered flux to that of the 
Rayleigh scattered flux 

RRam = FRam/ -^Ray 

(1 - rRay) p (3.10) 

rRay 

where 

Rs = y (3.11) 

For sufficiently large r^,, /3{n) decreases for small n. Because /(n) is decreasing suf- 
flciently fast, the ratio Rs is dominantly affected by the sum of the ffrst few terms. A 
conservative estimate can be obtained by considering only the first term, which gives 

_ (1 - rRay) 

RRam - ^(1) • (3.12) 
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Here, in the limiting case where Tg — > oo, can be evaluated exactly from the first 
principle, which is simply 

^ Wo io 8^ ^ ^ (3.13) 

= (11 - 121n2)/16 ~ 0.168. 

This is approximately equal to {^e)~^ given either by Eq. (3.5) or by Eq. (3.7). Hence for 

'^Ray = 0.8, which is close to the true value of the O VI doublets, we have -R^^^ = 1.48 
using This ratio is a typical one when the scattering optical depth is moderately 

larger than 1. 

For a slab of = 10, if we include terms up to n = 20, the ratio becomes RRam — 2.11 
which differs with -R^J^m significantly. Note that this number is sensitively dependent on 
the total scattering optical depth of the slab and that it usually increases with Tg. This 
implies that even though the Rayleigh scattering cross section is much larger than that 
of the Raman scattering cross section, the emergent Raman scattered flux is comparable 
to or larger than that of the emergent Rayleigh scattered flux, when the medium is not 
optically thin. This is also confirmed by Schmid (1992) and by Harries & Howarth (1997). 

For a small scattering optical thickness, the flux ratio of the reflected components 
is fairly close to the ratio of the scattering cross sections because the effect of multiple 
scattering is negligible. In the above formalism, this can be shown in a straightforward 
way. That is, f{n) ~ for n > 2, and we may take /3(1) ~ 1. Therefore, we have 

D d(1) _ 1 ~ fRay _ (^Ram /o a\ 

which implies that no conversion to Raman scattering from Rayleigh scattering occurs 
signiflcantly. 



4. Numerical Results 



In this section, we discuss the results of our Monte Carlo calculation for various cases 
including non-spherical wind models. The largest polarization is expected for 90 degree 
scattering, and negligible polarization is obtained for forward and backward scatterings. 
In this section the photons collected from the numerical simulations are those escaping 
perpendicular to the binary axis unless stated otherwise. Because of the symmetry about 
the binary axis in spherical wind models the polarization direction is either along the 
binary axis or perpendicular to it. Therefore we denote by negative P the polarization 
along the binary axis and by positive P the polarization perpendicular to it. 

4-1 Simple Gaussian Case 

The physical conditions of the emission line region are not known in detail, and we 
assume that it is in thermal equilibrium with a temperature T^mi being located very near 
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to the hot star. Under this assumption, the profile of the emission fine is described by a 
Gaussian with the width of 

vth = \^2kT^rn/ma = 2.3 Tg^^4 km s"\ (4.1) 

where Tgrn,A = T'em/(10^ K), and the oxygen atomic mass = 16 m^, rup being the 
proton mass. We take the wind terminal velocity v^o = 20 km s~^, = 50Rq and 
the separation of the cool giant and the hot star a = lOi?* as in Paper I. From now on, 
we will refer this case as the simple Gaussian case, which will serve as a reference model 
throughout the paper. 

In Fig. 6 we show the line profiles and the polarization for both A1032 and A1038 
photons for the cases where tq = 0.5, 1.0, 10. 

The polarization behavior and the line profiles are similar to those of the single scat- 
tering case, where the large polarization perpendicular to the binary axis is shown with 
small flux in the blue part of the feature and the polarization becomes small in the center 
and red parts of the feature. Furthermore, the polarization direction changes around the 
center part of the feature when Tq ^ 5. When tq ^ 1, the total flux shows a double-peaked 
proflle and in the case of opposite limit, a broad single-peaked profile is obtained. This 
is already summarized in Paper I and also by Harries & Howartli (1997) (see also Schmid 
1996). 

It is noted that the single scattering approximation gives a good qualitative description 
when To ^ 10. Another important point to note is that the scattering region corresponding 
to the polarization fiip is the region which is characterized by the Doppler factor about 
zero, that is, the sphere having the binary axis connecting the hot component and the cool 
giant as a diameter. We may give a simple intuitive explanation for this phenomenon in 
the context of the single scattering approximation as follows. 

The boundary region for the polarization direction fiip is the conical region with the 
opening angle of 7r/4 with the apex at the hot component. From Fig. 2, the scattering 
region responsible for the blue shift is mostly inside the conical region giving strong po- 
larization perpendicular to the binary axis. However, in the center-red part there is a 
competition of the two polarization components because the region with Doppler factors 
and larger lies both inside and outside the conical region. The contribution outside the 
conical region increases as the typical scattering optical depth tq becomes large. This 
explains the phenomenon that the polarization flip around the central region is obtained 
clearly when tq is signiflcantly larger than 1. 

4-2 Extended Ionized Region 

The assumption that the UV emission line region is spatially point-like is rather naive 
and the scattering region reduces severely due to the extended emission line region. Several 
researchers investigated the detailed shape of the emission nebula around the hot compo- 
nent of a symbiotic system. Seaquist et al. (1984) provided the shape of the ionization 
front using a parameter Xh deflned by 

^ ( 1^''] (4.2) 
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where Lh is the total luminosity of the ionizing radiation, M is the mass loss rate of 
the cool giant, v is the wind speed, a is the separation of the cool giant and the hot 
component, and as is the case B recombination coefficient for hydrogen. Monte Carlo 
calculations for various values of Xh have been done for the Rayleigh scattering of UV 
radiation (Schmid 1995) and also for the Raman scattering of the O VI lines (Schmid 
1996, Harries & Howarth 1997). We show some of our results to point out the profile and 
polarization features affected by the parameter Xh- 

We consider three cases where the parameter Xh is given by Xh — 0.04, 0.4, 4.0. 
For a value of Lh = ISOLq, where Lq is the solar luminosity, the corresponding mass 
loss rates of the cool giant are M = 10~^, 10~^, 10~^ Mq yT~^. In Fig. 7 is shown a 
typical ionization structure around the hot star. The solid lines represent the ionization 
front for the H II region, and the dashed lines are the conic sections which we adopt to 
approximate the true location of the front {i.e. the solid lines) in order to speed up the 
numerical calculation. 

The O VI line photons originate deep inside the H II region, because O VI has much 
higher ionization potential than H I does. Therefore, we assume that the O VI region is 
well- localized around the hot star, having a small size compared with that of a typical 
scattering geometry and justifying the point-source assumption. The results are presented 
in Fig. 8, where we only show the cases tq = 1, 10 for the A1032 photons. 

The overall results are qualitatively similar to those of the simple Gaussian case de- 
scribed in the preceding section. However, when Xh gets as large as ~ 4, the polarization 
is dominated by the component perpendicular to the binary axis, and no polarization flip 
is seen around the center of the feature. A simple explanation of this behavior is that the 
extended H II region has more intersection with the scattering region responsible for the 
parallel component of the polarization than the region giving the perpendicular compo- 
nent in the simple Gaussian case (c/. section 4.1). Therefore, compared with the simple 
Gaussian case, the reduced scattering region for the parallel component results in the per- 
pendicular polarization on the whole. This tendency is intensified as Xh increases. Hence 
combined with the analysis of section 4.1 the polarization perpendicular to the binary axis 
is stronger as Xh increases and tq decreases, which is displayed by the lower left panel of 
Fig. 7. 

4.3 Disk Wind 

As a non-spherical model a disk wind is invoked and described in this subsection. For 
ease of Monte Carlo coding, we just exclude the conic section of the spherical stellar wind 
from the scattering region. That is, the stellar wind is described by 

Vd..fc(p, 9y) = I - P~'^'^ 'I ^ ^ - (4.3) 

to if otherwise 

where 9o is the opening angle of the disk wind. Here, the normal direction of the disk 
is chosen to be the y-axis and the 2;-axis coincides with the binary axis. The angle 9y is 
defined as the angle making with the j/-axis, i.e. , 

9y = tan-1 ^^±Z. (4.4) 
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The density law is given in a similar way as in Eq. (4.3). 

In Fig. 9 we display the numerical results from the disk wind with an opening angle 
^0 = 7r/4. The results are displayed separately for the two cases. In the first case the 
observer sees the symbiotic system from the equatorial plane (Fig. 9a) and in the other 
case the observer is on the polar direction (Fig. 9b). Even though the azimuhtal symmetry 
about the binary axis is broken, the dominant polarization turns out to be either along the 
binary axis or perpendicular to it, and hence we keep the convention of the polarization 
sign to denote the direction of polarization. 

In both the cases the profiles and the degree of polarization are qualitatively similar to 
those of the simple Gaussian case. The plausible explanation for the results is that in the 
range where the single scattering approximation is a valid one the same Doppler-Tg contour 
diagram (Fig. 2) applies and therefore a similar argument leads to the qualitatively same 
results. 



As another application of non-spherical wind models, we investigate a bipolar wind, 
which is simplified as a truncated spherical wind model in a similar way to the disk wind 
case described in the previous subsection. The wind symmetry axis is assumed to be 
perpendicular to the binary axis. Therefore, the velocity law is 



The numerical results are displayed in Fig. 10, where the wind opening angle = 7r/4. 
In the left panel, the observer's line of sight is perpendicular both to the bipolar wind 
direction and to the binary axis, whereas in the right panel the observer lies in the direction 
of the wind. 

Unlike the disk wind case, we do not see any polarization flip in the bipolar wind 
model. This is again explained in a similar way to the case of the extended ionized region, 
where a large fraction of the scattering region responsible for polarization flip in Fig. 2 is 
excluded in the bipolar wind model. Therefore, the overall behavior of the polarization 
is characterized by a large degree of polarization in the blue part of the feature and the 
existence of residual polarization and flux in the red most part. 



Some of symbiotic systems exhibit complex features in their emission lines. The de- 
tailed physical nature of the emission line region of a symbiotic system is controversial and 
the velocity scale associated with the stellar wind around the hot star is sometimes of order 
10^ km s~^{e.g. Vogel & Nussbaumer 1994). Osterbrock(1970) proposed that synthetic 
line profiles composed of triangular shape may approximate the observed profiles quite well 
in the planetary nebula IC 418. We also use single-peaked and double-peaked triangular 
synthetic profiles to compute the profiles and polarization of the Raman scattered lines. 

The parameters necessary to describe the profiles are the widths and the peak values. 
In this work, the half width at the bottom of the single-peaked triangular profile is assumed 



4-4 Bipolar Flow 




(4.5) 



4-5 Synthetic Profiles 
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to be 30 km s~ . For the double-peaked triangular profile we superpose the two single- 
peaked triangular profiles with different peak values. The ratio of the peak values in the 
double-peaked profile is taken to be 0.3. 

In the right panel of Fig. 11 are displayed the Monte Carlo results for a single-peaked 
profile. The input profiles are shown in a small box at upper right corner in each panel of 
Fig. 11. As stated in Paper I synthetic profiles can be decomposed into 5-function profiles, 
which is the basis of the Green function formalism. In this regard, we see a large degree 
of polarization of the blue part of the feature, reminiscent of the simple Gaussian case 
discussed in section 4.1. 

In the left panel of Fig. 11 we show the corresponding results for a double-peaked 
profile. There are two peaks in the degree of polarization corresponding to the blue wing 
of each triangular peak. Therefore, the locations of the peaks in the degree of polarization 
shift to the blue with respect to those of the scattered fluxes. The amount of the shift 
depends mainly on the wind terminal velocity, and it increases as the wind velocity gets 
larger. Hence, the wind terminal velocity can be constrained by the shift of the peaks in 
the polarized fiux with respect to the total scattered fiux. 



It is expected that the hot component and the cool giant of a symbiotic star revolve 
around their center of mass with a velocity 



where ni — ///(I Mq) is the reduced mass of the hot component and the cool giant in 
terms of the solar mass and ai = a /{I AU) is the separation of the two components in 
the astronomical unit. This velocity scale is comparable to the terminal wind velocity 
and therefore may affect significantly on the profiles and the polarization of the Raman 
scattered features. However, it is known phenomenologically that the typical separation 
a of the 'S' type symbiotic stars is very different from that of the 'D' type symbiotics 
{e.g. Iben Sz Tutukov 1996). In this respect, the spectropolarimetry monitored for an 
extended term will reveal a good deal of information about kinematics. In this section, we 
investigate the effect of the orbital motion of the hot component with respect to the cool 
giant on the profiles and the polarization of the Raman scattered features. 

In Fig. 12 we draw the Doppler factor-Tg contours analogous to Fig. 2, where the hot 
component is moving in —y direction with velocity Vrev = 0.6 Voo, and the observer's line 
of sight coincides with x-direction. The Tg contour remains the same because it is assumed 
that the wind structure around the cool giant remains spherically symmetric in the rest 
frame of the cool giant and the scattering cross section is insensitive to the velocity scale of 
order Vqo = 20 km s~^. The diagram is plotted in the rest frame of the cool giant. Due to 
the enhancement of the Doppler factor, the relative motion of the cool giant with respect 
to the observer is not important compared with the motion of the photon source relative 
to the scatterers. 



5. Effect of Binary Motion 




(4.6) 
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Notable changes are seen in the structure of the Doppler factors. The contour corre- 
sponding to zero Doppler factor is not spherical any more and the sphere connecting the 
two stars may be decomposed into the upper hemisphere with positive Doppler factors and 
the lower hemisphere with negative Doppler factors. The important point to note about 
this sphere is that in the case of simple Gaussian case (see section 4.1) the polarization flip 
occurs around this sphere, which corresponds to the central part of the scattered feature. 
Therefore, we can expect from this diagram that there is possibility of polarization flip in 
the red part and/or the blue part of the scattered feature. 

The numerical results corresponding to the situation depicted in Fig. 12 are presented 
in Fig. 13, where the observer's line of sight is perpendicular to the orbital plane. Here, we 
also present the position angle because the symmetry about the binary axis is broken with 
the introduction of the relative motion of the hot component with respect to the cool giant. 
The position angle is not defined when the degree of polarization is 0. Therefore in the 
extreme red part where there is so small an amount of flux not to be assigned measurable 
polarization, the position angle is highly uncertain and does not yield any quantitative 
information. 

When the scattering optical depth tq = 1 and v^ev = 0.6 Voo, the polarization flip 
occurs around the red part. The degree of polarization is strong in the blue part of the 
feature. The position angle remains almost perpendicular to the binary axis in this part 
and parallel in the red part of the feature. The polarization component perpendicular to 
the binary axis forms a double-peaked structure and therefore with the parallel component 
the overall polarization shows a triple-peaked structure. The triple-peaked structure is also 
seen in the polarized fliix. 

In the total scattered flux no obvious multiple structures are seen in these cases. The 
spectropolarimetry provided by Harries & Howarth (1996) shows that many symbiotic 
stars also have the multiple peak structures in the flux as well as in the polarization and 
the polarized flux. In the simple Gaussian case with small optical depth (tq < 0.5) there 
is a faint double peak structure in the fliix. However, the small parameter space is not 
enough to explain the much more complicated phenomena shown by the observation. This 
strongly implies that the simple Gaussian proflle of the incident radiation is difficult to 
explain the complicated behavior of the spectropolarimetry. 

It is also noted that the polarization flip occurs at a scattering optical depth tq = 1. 
This is in contrast with the simple Gaussian case, where the polarization flip at the center 
part of the feature is obvious only when tq ^ 5. This is particularly interesting consid- 
ering the spectropolarimetry data on the symbiotic systems including RR Tel, where the 
polarization flip is clearly seen in the red part of the Raman-scattered features (e.g. Espey 
et al. 1995, Schmid & Schild 1990, Schild & Schmid 1996, Harries & Howarth 1996). 

When the orbital velocity exceeds the wind terminal velocity, then the polarization flip 
is also shown in the blue part, which is illustrated in Fig. 14. Here, we adopt Vrev = 3 Vqo- 
The polarization flip is seen in the red part and also in the blue part with small flux. 
This behavior may be explained using a similar argument mentioned in the preceding 
paragraphs with Fig. 11. 

In Fig. 15, we show the numerical results when the observer's line of sight lies on the 
orbital plane. The results are similar to those of the simple Gaussian case, and therefore 
no polarization flip around the red part occurs in this case. This dramatically different 
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behavior can be also explained using the Doppler factor-Tg diagram. When the observer's 
line of sight lies on the orbital plane, the section of the Doppler factor- diagram including 
the binary axis and perpendicular to the observer's line of sight coincides with Fig. 2, and 
the polarization behavior is not significantly different from that of the simple Gaussian 
case. 

6. Observational Implications 

Observational and evolutionary properties of symbiotic stars are extensively summa- 
rized by Iben & Tutukov (1996), according to whom symbiotic stars are rather inhomo- 
geneous group showing a large range of variabilities and physical conditions. Detailed 
analysis of the UV emission lines of symbiotic stars also implies that kinematics around 
the emission line region may not be described by a simple unique model. 

In the optically thin limit, the ratio of the A1032 flux to that of the A1038 flux in 
the O VI doublet is 2, which is the ratio of the oscillator strengths of the two transitions. 
However, this ratio is altered when the medium gets optically thick, where the stronger 
lines have more difficulty escaping the region. Far UV observations show that for the 
doublets including C IV AA 1548, 1551, N V AA 1238, 1241 the ratios are usually less 
than 2, and sometimes less than 1 {e.g. Vogel & Nussbaumer 1994), which conflrms the 
complicated nature of the emission line region. 

The Raman scattered features provide a special tool for investigation of the physical 
conditions of a symbiotic star through spectropolarimetry. Schmid & Schild (1994) pro- 
vided spectropolarimetric data of some symbiotic systems, and more recently Harries & 
Howarth (1996) presented more data with enhanced resolution. The polarization behavior 
is also very heterogeneous and has complex structures. Schmid & Schild (1994) divided 
them into the three types and described the basic points. The extended ionized region 
may have definite effects on the polarization behavior displayed in the observational data 
as shown by Schmid (1996). 

It has been pointed out by several researchers including Schmid (1996), and Harries 
& Howarth (1996,1997) that the position angle varies throughout the Raman scattered 
features. This implies that the most general model should include those producing the 
polarization direction other than the binary axis. Because the polarization behavior based 
on the models possessing a symmetry with respect to the binary axis does not show any 
rotation of the position angle, it has been proposed to consider non-spherical models. 

Harries & Howarth (1996) pointed out that most of the Raman scattered features 
show multiple peak structures in the flux and in the polarized flux, and in particular 
double and triple peaks are found very frequently. These structures are maybe coupled 
with the complicated proflles of the UV emission lines including a P-Cygni type proflle. If 
the incident UV line proflle has a double-peaked structure, the similar structure is obtained 
in the scattered proflle and in the polarized flux. However, the location of the peaks may 
differ depending sensitively on the wind terminal velocity. In this regard, it is necessary to 
consider various input proflles for reproducing the characteristics of the observed Raman 
scattered fluxes. 

Some fraction of symbiotic stars show polarization flip around the red wing of the 
scattered features, which is produced in the Monte Carlo computation when the orbital 
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motion of the hot component about the cool giant is included. Typical parameters of 
symbiotic stars imply that the rotation velocity is almost comparable to the wind terminal 
velocity {e.g. Eq. (4.6)), and the polarization flip is shown clearly when the orbital plane 
coincides with the sky plane. The orbital motion can be observed by a careful monitoring 
of the position angle variation {e.g. Schild Sz Schmid 1996). However, we predict that the 
polarization flip will be seen irrelevant of the orbital phase, because the inclination angle of 
the orbital plane with respect to the line of sight does not change with the orbital motion. 

Furthermore, the polarization flip is shown in the blue part when the orbital velocity 
exceeds the wind terminal velocity. In this case, the polarized flux shows a clear triple- 
peaked structure. Therefore, an independent observational determination of the orbital 
velocity will provide a good constraint to the wind terminal velocity, which is proportional 
to the mass loss rate of the cool giant. 

It is proposed that the effect of the binary motion may leave an important signature in 
the polarization and polarized flux. It will be also interesting to investigate the combined 
effects of the binary motion and the non-Gaussian input profiles with different ionization 
structures. However, the introduction of the binary motion breaks the azimuthal symmetry 
and requires a good amovmt of computation time and a huge volume of parameter space 
to be examined. Therefore, independent determination of the kinematic and the dynamic 
parameters by the UV observations and monitoring of spectropolarimetry are expected to 
unveil the physical nature of symbiotic stars. 

7. Summary 

The Raman scattered features are quite unique to symbiotic stars so far and can be 
used as a powerful diagnostic of the physical conditions including the ionization structures 
and the kinematics on the orbital motion and the stellar wind structures. They possess 
several distinguished characteristics. The first and the most important point is that they 
are usually strongly polarized and therefore spectropolarimetry reveals a significant amount 
of information about the kinematics and the scattering geometry. Related with this is 
that they are composed of purely scattered photons without dilution from the direct flux 
retaining the vivid information on the scattering geometry. The third point is that the 
inelastic nature of the Raman process provides the enhancement of the profile width by a 
factor of A//Ai, which amounts to almost an order of magnitude in the case of the O VI 
doublets. 

In this study, we reproduce the general properties of the Raman scattered emission 
lines, which are found in the works of many researchers including Schmid (1995, 1996) 
and Harries & Howarth (1997). The spherical wind model usually gives a large degree of 
polarization in the blue wing and the polarization fiip around the center of the feature 
depending on the scattering optical depth. 

If we take the orbital motion of the hot component around the cool giant into consid- 
eration, then the polarization fiip is seen to shift to the red wing when the orbital plane is 
near to the sky plane. This effect is conspicuously seen when the orbital velocity is almost 
comparable to the terminal velocity of the stellar wind, giving insight into the spectropo- 
larimetric data of several symbiotic stars including RR Tel {e.g. Espey et al. 1995, Harries 
& Howarth 1996, Schmid 1996). When the orbital velocity exceeds greatly the terminal 
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velocity, the polarization flip is expected to occur both in the blue wing and in the red 
wing. We also expect that various input profiles proposed by generic dynamical models 
can be combined to produce the complex features in the Raman scattered lines. 

Much parameter space remains to be investigated by numerical methods, including 
the parameters M, Xh, Voo, Vrev to list a few, and it is hoped that the numerical 
calculation may provide a useful constraint to existing dynamical models. It is generally 
concluded that spectropolarimetry will remain a useful tool to delve into the symbiotic 
system. 
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FIGURE CAPTION 



Figure 1.- A schematic geometry of a symbiotic star system for the Monte Carlo 
calculation. The 2;-axis is chosen to be the binary axis and the observer's line of sight 
coincides with x-axis. The hot star is represented by the circle on the right side, and 
the bigger circle in the origin is the cool giant. A spherically symmetric stellar wind is 
shown by radial (small) arrows. The vectors pi and pi+i represent the i-th and the i + 1-th 
scattering sites respectively and siji) is the spatial distance between the two scattering 
sites, h is the impact parameter of the photon trajectory with respect to the cool giant. 

Figure 2.- A contour map of the Doppler factors and the scattering optical depths 
of incident photons from the hot star. The stellar wind is a spherical one defined in 
section 2.2. The Doppler factor is defined by DF = kj • v/c and the total scattering 
optical depth = j ds n{s){aRay + (^Ram) ■ The adopted parameters are described in 
Paper I. 

Figure 3.- Comparison of the single scattering approximation with corresponding 
Monte Carlo calculation. The dotted lines represent the Monte Carlo result and the result 
from the single scattering approximation is shown by the solid lines. 

Figure 4.- The Rayleigh and Raman reflected fluxes from a finite slab of total scattering 
optical depth Ts = 0.5, 1, 5, 10 as a function of the number of scattering. The light lines 
represent the result for the case RRay = 0.2 and the case for Rray = 0.8 is shown by the 
thick lines. The solid lines are for the Raman scattered fiiix and the dotted lines are for 
the Rayleigh case. 

Figure 5.- The Rayleigh and the Raman reflected fluxes from a slab of the total 
scattering optical depth Ts = 10, RRay = 0.8 as in Fig. 4(d) with semi-quantitatively 
described fRamin), fRayin) in section 3.2. The solid lines are for the Monte Carlo results, 
and the dotted lines give fRamin), fRayin) computed using Eqs. (3.3)-(3.6). The upper 
set of curves are Raman scattered lines and the lower set corresponds to the Rayleigh flux. 
In the bottom panel, fRamin), and fRay{n) calculated from a different escape probability 
(3'{n) = exp(-n0-4)/2 (Eq. (3.7)). 

Figure 6.- The line proflles and the polarization for the simple Gaussian case, where 
the UV emission line proflles are described by a Gaussian with Tgm = 10^ K. See the text 
for more detail. 

Figure 7.- A typical ionization structure around the hot star of a symbiotic system, 
reproduced using the recipes provided by Seaquist et al. (1984). The parameter Xh deflned 
in Eq. 4.2 determines the overall structure of the ionized region. Solid lines divide H I region 
for H II. Dotted lines are obtained using a simpler equation for the ionization front. 

Figure 8.- Monte Carlo results corresponding to the ionization structures described in 
Fig. 7. 

Figure 9.- Monte Carlo results for the disk- wind models. 
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(a) The observer's line of sight hes in the equatorial plane of the disk wind and the 
typical scattering optical depths ro = 0.5, 1.0, 10. 

(b) The observer's line of sight is perpendicular to the equatorial plane. 

Figure 10.- Monte Carlo results for the bipolar wind models. 

(a) The observer's line of sight is perpendicular both to the binary axis and to the 
wind symmetry axis and the scattering optical depths tq = 0.5, 1, 10. 

(b) The observer's line of sight coincides with the wind symmetry axis. 

Figure 11.- Monte Carlo results for the cases where the input profiles are gives by 
single-peaked and doublc-pcakcd profiles. The input profiles are shown in a small box at 
upper right corner in each panel. 

(a) Single-peaked profile with tq = 1 

(b) Single-peaked profile with tq = 10 

(c) Double-peaked profile with tq = 1 
(b) Double-peaked profile with tq = 10 

Figure 12.- Doppler-r^ contours analogous to Fig. 2 including the relative motion of 
the hot component with respect to the cool giant. The hot star is assumed to move in — 
direction with a velocity Vrev = 0.6 v^o, where Vqo is the terminal velocity of the spherical 
stellar wind around the cool giant. 

Figure 13.- Monte Carlo results corresponding to the case described in Fig. 12. The 
observer's line of sight is perpendicular to the orbital plane. 

(a) The typical scattering optical depth tq = 1 

(b) The typical scattering optical depth tq — 5 

Figure 14.- Monte Carlo results for with Vrev = 3 Voo ■ The wind structure around the 
cool giant and the observer's line of sight are the same as in Fig. 12. 
(a) To = 1, (b) To = 5 

Figure 15.- Monte Carlo results for the same geometry with Fig. 12, where the ob- 
server's line of sight lies on the orbital plane, 
(a) To = 1, (b) To = 5 
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